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Abstract 



An algorithm to construct analytic approximations to two-loop diagrams de- 
scribing their behaviour at small non-zero thresholds is discussed. For some 
special cases (involving two different-scale mass parameters), several terms of 
the expansion are obtained. 
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1 INTRODUCTION 



The calculation of two-loop Feynman diagrams with (different) masses meets serious technical difficulties, 
even for two-point functions. After the tensor decomposition of two-loop self-energy diagrams [Q, the 
problem is reduced to scalar integrals. When no exact analytic expressions are available, two main 
strategies can be chosen: either numerical or approximative analytical. As to the numerical strategy, there 
exist various integral representations |Q which provide results for given masses via numerical integration. 
A semi-numerical approach was considered in Q . 

The analytical approach uses explicit formulae for the asymptotic expansion of Feynman diagrams 
in momenta and masses and is based on general theorems on asymptotic expansions Q (see also ||^ 
for review). For two- loop self-energy diagrams with general masses, the first results have been obtained 
outside the particle thresholds: a small momentum expansion below the lowest threshold |q and a large 
momentum expansion above the highest threshold (some special cases were considered in^]). However, 
small and large momentum expansions do not describe the behaviour between the lowest and the highest 
physical thresholds. In the special case of a zero mass threshold the small momentum expansion can be 
extended to the lowest non- vanishing threshold |^ (see also in ||l^). The expansion coefficients involve 
then powers of ln(— A:^) (where k is the external momentum). 

Our investigation considers cases when one (or two) two-particle threshold(s) is (are) small with 
respect to the other thresholds, but not anymore zero. Using asymptotic expansions in the large mass 
limit, a series converging above the small threshold can be found. The expansion coefficients now contain 
the two-particle cut(s) associated with the small threshold(s) and therefore the non-regular behaviour 
around the threshold(s) is described. Here we restrict the applications to diagrams involving one large 
(M) and one small (m) mass parameter, i.e. some masses are equal. In | ]TT| a more detailed description 
of the approach and of the expressions for general masses is given. 



2 THE APPROACH 

Consider the diagram shown in Fig. 1. The corresponding scalar Feynman integral reads 

J({.J;W;fc) = | J (1) 

where n = 4 — 2e is the space-time dimension and [DiY^ = {p^ — mf -\- iO)''' are the powers of the 
propagators corresponding to the internal lines. The momenta pi are constructed from k and the loop 
momenta. 

In general, the diagram in Fig. 1 has two two-particle thresholds, at fc^ = (7711-1-7714)^ and k^ = 
{m2 + 'rn^Y , and two three-particle thresholds, at = (rni + m^ + m^Y' a.Tid k'^ — (rn2 + '013 + 1714)'^. Let us 
consider some of the masses to be large, while the other masses and the external momentum are small. We 
shall denote the large masses with capital letters. Mi. The classification of all different small-threshold 
configurations has been given in refs. |ll| (see also in Table 1). 

We now need to introduce some notation. Let F be the original graph (Fig. 1), its subgraphs are 
denoted as 7, and the corresponding reduced graph r/7 is obtained from F by shrinking 7 to a point. 
Furthermore, is the dimensionally-regularized Feynman integral with the denominators corresponding 




Figure 1: The master two-loop two-point diagram 
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Table 1: Subgraphs 7 to be included in the expansion. 



Case 


Heavy 
lines 


Subgraphs 7 


1 


1, 3, 4 


r,{i34} 


la 


1,4 


r,{1245}, {134}, {14} 


lb 


1, 3 


r,{1235}, {134}, {13} 


2 


3 


r,{134}, {235}, {3} 


3 


1, 5 


r, {1235},{1345},{1245},{15} 


4 


1 


r,{1245}, {134}, {1} 



to a graph 7. In particular, Jr corresponds to (|^) itself. Then the general theorem [Q yields 

Jr '^Jr/l°1k,nH,qiJf, (2) 

fc, TTli^O ^ 

where the sum goes over all subgraphs 7 which (i) contain all the lines with the large masses Mi, and 
(ii) are one-particle irreducible with respect to the light lines (with the masses rrii). The Taylor operator 

mi, qi expands the integrand of J-y in small masses m^, the external momentum k and the loop momenta 
Qi which are "external" for a given subgraph 7. The symbol "o" means that the polynomial in qi, which 
appears as a result of applying T to J-y, should be inserted in the numerator of the integrand of Jrf-y- 

Let us consider which subgraphs 7 contribute to the sum (||) for the different small-threshold config- 
urations, using the numbering in Fig.l. For example, the "full" graph F = {12345} includes all five lines, 
{134} corresponds to a subgraph without the lines 2 and 5, etc. The 7's contributing to the expansion 
are shown in Table 1. 

After partial fractioning, the following types of contributions arise in the coefficients of the small- 
threshold expansion (in situations with two-particle small thresholds): (a) two- loop vacuum diagrams 
with two (or one) large-mass lines and one (or two) massless lines; (b) products of a one-loop massive 
diagram (with small masses and external momentum k) and a one-loop massive tadpole; (c) products of 
two one- loop massive diagrams with external momentum k. For the diagrams containing a small three- 
particle threshold, we also need results for the diagrams obtained by shrinking the lines corresponding to 
large masses. 

The contributions of type (a) are discussed in 1^. The one- loop two-point functions involving small 
masses of internal lines (occurring in the contributions (b) and (c)) are "responsible" for describing the 
two-particle threshold irregularities (some relevant results are collected in [pT|). 

3 SOME RESULTS 

As an example, consider the diagram in Fig. 1 with 1^1= ... =vq = 1 and n = 4. Denote the corresponding 
integral (|I]) as (here, rui may correspond to either small or large masses) 

OQ 

J(mi, TO2, ma, 7714, ms; fc) = — tt'' ^ S'j , (3) 

j=o 

where Sj are the terms of our expansion. For a given j, the term Sj is a sum of all contributions of the 
order (fc^)^° IK™?)"'' (where the product is taken over all lines with small masses) with jo — J- 

The Sj may contain logarithms of the ratios of the masses and some functions which are described below. 

The contributions of type (a) involve a function 7i(M^, M|) which can be expressed in terms of Li2 
(see in 1^). It is antisymmetric and therefore vanishes when Mf = M|. 

The contributions of types (b) and (c) are expressed through the function r(7n,i, 7712; fc^) which is 
related to the finite part of the one-loop self energy and can be expressed in terms of elementary functions 
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Table 2: Terms of the expansion for the cases 1 and la 



Case 1: J(M, m, M, M, m; k) 

M'^So — + — I 

M^S, = e (-^r + ^L™ -^)+rn^ {-^^r + iL„ - §) 

M^Sq, = (P)^ ("TM''' + TM-^™ ~ Tm) + ^^^^ ("M''' + M^™ ^ 215) + (^M''" + ^-^"1 ~ ^) 
-^^^'3 = (P)^ ("Ti^''' + TT20-^™ ~ 448o) + {k^Yrm? {--^t + -^L^ - 19^) 

,l2 4/'_J_ I J_r _ ^9_\ I 6 J__ I J_r 743 \ 

^fi, /(t / gg-J^m 1568/^'" V 280 ' ^ la^rn 29400/ 

M^°S'4 = (^eM)''' + 6M)^™ ~ loSoo) + {k'^Yrm? {-^^t + g|o^m + 99525) 

_L('7,2\2 4 ^T-iJ_r 3067 \ I 7 2^6 I i J_ r 1333 \ 

-ri^A, y //t 1^ ^gpp / -I- jgg -t^m 1058400/ I, 1260 ' 126 -^"^ 158760/ 

+to8C ^T-i J_r _ 2131 \ 

1260 ' ^ 252^™ 317520/ 

Afl2q^ _ (i„2\b ( 1 ^1 1 r 1 \ I (].2\4 2 ( 1 ^1 I 1091 \ 

ivi 05 — ) y 33264' ^ 33264^™ 66528/ ^ v'*' / '"■ V 16632' ^ 2772"^™^ 7683984/ 

+ 4 / 1^1 53 r 6047 \ , (p)2 6 ( ^^-i ^3_r 191071 \ 

-ryi^ ) III. \ 9240 ' ^ 27720 ™ 48024900/ ^ > '"' \ 5940 ' ^ 5940 " 



27720 ™ 48024900/ ' / V 5940' ' 5940 ™ 82328400/ 

I 1,2^8 ( _i 5_r _ J2153_\ , 10 ( 1 J_ r 22727 \ 

-Th. in y ^^g2 ' "I" 1584-^™ 627264/"'"'"' V 5544' ^ 924-^™ 12806640/ 
Afl4c _ (U2\(, ( 1 -^1 1 r 1 \ I /l2\5 2 ( 1 -^1 1 r I 126421 \ 

— ; 168168' ^ 168168 ™ 392392/^'^'^ / '" V 84084' ^ 12012-^^^ 2164322160/ 

I /l2'|4^4 ( 1 ^ I 43 r 1 433253 1 /p')3^6 / 1 -^1 313 r 6920149 \ 

-r^n, y //I ^^ 45864' """ 72072^™""" 2597186592/ > '"" V 28028' ^ 180180^™ 16232416200/ 

I (1^2\2^?, ( 1 _ I 1 r 63169 \ , )U2^10 1^ 1 ^ , 1 r 100453 \ 

-ryt^ / V 20020' """ 468 »" 42162120/"'"'^ V 18018' "^ 858 77297220/ 

+rn}'^( ^—T+-^L 288851 \ 

\ 24024' ^ 3432-^"* 618377760/ 

Case la: J{M, m, m, M, m; k) 

M'^So = -T + L„ - 2 

M^Si = e {-^T + ^Lrn + T2) + m2 (-tL„ - |r + + 2L„ + 2C2 - 3) 
M'S, = {kr (-^r + + ^) + fcW (ir + L„ + f) 

+7X1" {-ArLm - f r + + 12i„ + I2C2 - f ) 

+Pm4(iTL„ + fr+|L^ + f L,„ + 3C2 + f )+m6(-15TL„-fT+30L^+64L„+60C2-ip) 



Miss's 



— I \ 3150 ' ^ 3150-^™ ^ 63000/ ^ '"■ V 1050 ' "^ 42 "^ 4200 j 

+ (fc2)W(-2^r + 2L2„ + fL„+4C2 + ^)+fc2rn6(4ri„ + ifr+20L^ + 110L„+40C2 + 2|f) 
+m« (-56tL„ - iif r + 140L^ + ^-fL^ + 28OC2 - 

_ (•p^5 f ^T--l I 2531 \ I f/.2'|4 2 / 1 , 209 r i 39887 \ 

! \ 16632 ' ^ 16632-^™ ^ 997920/ ^ > '"' V 7560 ' ^ 1260-^™ ^ 75600,' 

(fc^)'™' (-io^ + + + 4C2 + S) 



+ , ■■- y 42y • ■ m • 140 — ' ' 5U40 / 

I /'^.2^2_,6 (' l_-r 29 , 115r2 , 3649 r , 230 ^ , 4989 \ 

+/c2to8 (f ri„ + i|^T + 175L^ + + 350(2 + ^) 

+mio (-210tL„ - if^r + 630i^ + 1459L„ + I26OC2 - ^) 

/\fl4o _ /; 2^6 <^ 1 ^1 ^ T 1 4903 \ , /; 2\5^2 1 , 113 r , 395999 \ 
ivi — \l^ ) y g^Qg^ / -I- 84084^™ ^ 3923920/ ^ ^-'^ / V 48510 ' ^ 990 ™ ^ 970200/ 



~ ) \ 84084' "•" 84084^™ """ 3923920/ "*" ' '"' V 48 

+ikrm' (- 4^10- + + W^"^ + 4C2 + ^) 

+ (fc2)3^6 + 562.2^ + l^i^ + 112C2 + 

+ (-4rL™ - fir + 462L,2„ + 12^^^ ^ 934(2 + 

+A:2mio (120tL„ + + 1260L^ + 5180i„ + 2520(2 + ^) 

+m^^ (-792rL^ - + 2772L^ + + 5544(2 - '-^] 
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Table 3: Terms of the expansion for the case lb 



Case lb : J{M, m, M, to, to; k) 

A'PSo = -T + L„, + C2 - 4 

M^Si = fc2 1^ + 1^^^^ + " if ) + - |r + + |i„, + 3C2 - ^) 

M6^2 = (fc')' (-|r + ii„ + iC2 -§§0)+ (-^^™ - + + li^m + 4(2 

+to4 (-3rL,„ - + 4L^ + + 9^ - W) 

M^S, = (fc2)3 (_^^ + iC2- Mi) + {k^m^ {-tL^- ffr 



12973 ') 
1800 / 



-k'^in^ 



47, 



15 



_ r2 

2 m 



I 2 6 

+to6 (-8tL„ - + 14L2„ 

= (^')M-2^-+^^™ + K2- 



5953 r I 39 A 
490 -'^m ~r o 1,2 



5953 r 
' 420 

1535 T I 00/- 73547 \ 



120 
1549283'! 
58800 / 



183 T 
70 ' 



5C2- 



765467 ' 
88200 . 



84 
86939 
252000; 



1764 / 



75 ' -"m 5Q4 

-TO« (-22tL„ - + 49L2 . 39869 

(fc2^5 / 1 



_ 32459123 '! 



+fc2TO6 (-20rL, 



504 



3175200 
o(;a _ 17514967'! 
"-""''2 396900 / 

QAA _ 138737377'! 

1260 ' ^"^"^ 1587600 / 

+ 99C2 



36 



Lt 4- I/- — 2828267 '! 
^fi-^m f;t,2 9979200/ 
, 49103 , 



41100499 '! 
317520 / 



36 



+ (fc2)3TO4 ( 



27T-r 

/' 136_r 
\ 



+k^m' 

+to1° (-64ri™ - 

+ (fc2)5„l2 ^ 

+(^2)3^6 (-92rL^ 
+ (^2)2^8 (-176rL 
+fc2TOio ('-192rL,„ 



630 

_ 16409 
560 
_ 62402, 
' 945 ' 
1069639,, 



13860 ^™ + '''^2 
35^2 I 18179 £ 



15120 
1943093, 
27720 



2 m 
I 283 7-2 

206Li 



360 
2806913 



175L2 



• ^ 10395 
I 9873841 

18480 
8884391 T 



4534026869 '! 
384199200 / 
ROA 20105059'! 
0U<'2 - 285120 j 

I 999/- _ 21361520641 ' 
™ "t" ^^^1,2 144074700 , 



144074700 
r 4- 849/' _ 6628516487 
2 ''2 21344400 



^ j_ 1 T _J_ I/' 85046849 



27720 



35IC2 



81848153249 ' 
192099600 / 



7689 , 
1960 ^ 
_ 829189, 
17640 ' 
557275, 



353152800, 
49963 



12870^™ 

24L 



■8C2 



3528 
6228527 



2955811 . 



13860 

12 (-196Ti„, - ^ffii^T + 637LI 



435006955477'! 
32464832400 / 

I qoA _ 7029366464489'! 

360360 ^™ ~'~ ^2 64929664800 / 
lSfir2 I 221210659 t i 479/- 14327879597579 '! 
lOUiv^ -I- 3gQ35Q 64929664800 / 

I -1 07c; A 5163216185287 '! 
-t-iO(0(,2 12985932960 / 
innfir2 , 4623547 r lOmAA 3647519729639 \ 
iUU»7v„, + 1^,38 ^rri + ^Ud4(,2 2951348400 ) 

L,„ + 1275(2 6314510897209 - 



2 _|_ 27760549 



__i_cc;nr2 I 756802831 t 

27720 ' "'^"^m 360360 



1638 
6324508 



5005 



4328644320 



(see e.g. in ||T^). By using the notation 

A = Am\m\ - {k^ 

we get 



TO 



2n2 



(4) 



t(toi, TO2; fc2) 




fc2_. 



+ (to2 — TO2) In ■ 



-i7rV^e'(fc2 - (toi +TO2)2)|. 



(5) 



111 our contributions, the t function depends on the small masses only. In fact, it contains the main 
information about the small-threshold behaviour at two-particle thresholds. The 9 term in the braces 
yields an imaginary part in the region beyond the physical threshold (i.e. for k"^ > (toi -|-TO2)2). Further 
properties of the r function are discussed in ||Tl|] . 

In Tables 2, 3, 4 we present results for the terms of the small-threshold expansion (||) for the cases 1, 
la, lb and 2 (cf. Table 1), provided that all large masses are equal to M and all small masses are equal 
to TO. We use the notation 



1^2 

gTT . 
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Table 4: Terms of the expansion for the case 2 
Case 2 : J(m, m. M, m, m\ k) 

A'PSo = -t' + 2tL,„ - 2t - L'^^ + 2L,„ - 2C2 

M^^i = fc2 (ir2 - ri,„ + ir + ii^^ - iL„, + (2 - |) + (-2t2 + 6rL,„ - r - 4L^ - 8(2 + lO) 
M^S^ = {ef {-\t^ + f rL™ - |r - + |L„ - §(2 + l) 

+Pm2 (f - f rL„ + fr + - 2L„ + 8C2 - W) 

+^4 (_^r2 + 20rL„ + f r - 18L^ - 14L„ - 86(2 + 47) 

M^S, = (fc2)3 (1^2 _ + 1^ + 1 ^2^ _ 1 + 1 _ 115) 

+ (fc2)2TO2 (_3^2 ^ 7^^^ _ 3^ _ 4^2^ 2^^ _ ^ ^) 

+fc2TO4 (12^2 _ 35rL„ - + 22^2^ - f L„ + 44(2 - ^) 

+^6 (_i6r2 + 70rL„ + if^r - 80L2„ - 292^^ _ 16OC2 + ^f^) 

= (P)4 (-ir2 + |rL„ - - + _ 2^2 + §) 

+ (fc2)3^2 (^^2 _ m^L^ + 16^ + 4^2^ „ ^^^^^ + _ 2^) 

+ (fc2)2^4 (_^^2 _^ 252^^^ ^ 58^ _ 34^2^ _ _ ^^^^ ^ 3187) 

+ P™6 (256^2 _ i68rL„ - + lOO^rn " + 200C2 " ^) 

+^8 (_256^2 _^ 252tL„ + - 350L2„ _ ifl^^ _ 700C2 + ^) 

Afl25'j. ^ (^2-)5 ^1^2 _ 1^-2,^ + + \L'^m " ll'^'m + K2 ~ iMo) 

+ (A:2)W (_^^2 + ^^^^ _ 1^ _ 4^2^ + - 8C2 + ^) 

+ (^2)3^4 (^^2 _ gg^^^^^ _ 2r + 38L^ - M2L,„ + 76C2 - 

+ (fc2)2^6 (. 320^2 _^ 308^L„ + 440^ _ 712^2^ _ 4474^^ _ 1^4^^ + 72|I^) 

+^2^8 (^r2 - 770tL„ - + 385L^ - if^L™ + 770C2 - ^) 

+^10 (.512^2 ^ 924^i^ + - 1512L2„ _ }m±L„^ _ 3024(2 + ^) 

/i/rl4c _ /^,2^6 / l_.2,2 r 2 lr2,2r 2^,139^ 

+ (fc2)5^2 (24^2 _ ^^^^ + 24^ 4^2^ _ 2^^^ 8C2 " ^) 
+ (^2)4^4 (_242^2 + 572^^^ ^ _ 50^2^ _ 2J2^^ _ ^qq^^ + 1^) 

+ (fc2)3^6 ( i^^2 _ ZA^^L^ - + 256Li - + 5I2C2 - iSfp) 

+ (fc2)2^8 (_3M0^2 + i7i6^2.,„ + - 1862L^ - - 3724(2 + ^) 

+kV (6M4^2 _ 3432^L„^ _ + lllQLl - I^L„, + 2352(2 - 2™) 

+^12 (_4cp^2 ^ 3432^L,„ + 22|p^ _ g4g8^2^ _ i^L„ - 12936(2 + ™23) 



In the Hmit m ^ 0, r {Lm — L), where L = ln(— fc2/M2). In this hmit the resuhs presented in 
Tables 2-4 reproduce the first four columns of Table 1 from ref. j^. To perform the calculations, we used 
the REDUCE system @. 

4 CONCLUSIONS 

We have studied the behaviour of two-loop self-energy diagrams when the external momentum and some 
of the masses are small with respect to the large masses. By use of explicit formulae for the terms of 
asymptotic expansions in the large mass limit, we presented an analytic approach to calculating these 
diagrams by keeping the first few terms of the expansion. The main idea was to avoid putting any 
conditions on relative values of the external momentum squared and small masses. In addition to the 
general results presented in here we listed the contributions (from 5*0 to 5*6) for some special cases. 
Some of the diagrams considered are interesting from the physical point of view, since they may occur in 
the Standard Model calculations. 

In cases 3 and 4, three-particle small thresholds arise. The small-threshold behaviour is then defined 
by the functions corresponding to the sunset diagram (with three propagators) and the diagram with 
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four propagators (e.g. eq. (|^) with i^i =0). Unfortunately, sufficient analytic information about these 
diagrams is not yet available. 

The algorithm can be also extended to the three-point two-loop diagrams with small thresholds (for 
zero- mass thresholds, see ref. jfSl). 

Acknowledgements. A. D. is grateful to the organizers (DESY-Zeuthen) for support of his par- 
ticipation in the Rheinsberg workshop. A. D. and V. S. were partly supported by the RFBR grant 
96-01-00654. 
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